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Abstract 



We derive source terms for the production of quarks and gluons from the 
QCD vacuum in the presence of a space-time dependent external chromofield 
Aci to the order of S^^\ We found that the source terms for the parton pro- 
duction processes Ad qq and Ad, AdAd gg also include the annihilation 
processes qq Ad and gg Ad, AdAd- The source terms we derive are 
applicable for the description of the production of partons with momentum 
p larger rhan gA which itself must be larger than Aqcd- We observe that 
these source terms for the production of partons from a space-time depen- 
dent chromofield can be used to study the production and equilibration of 
the quark-gluon plasma during the very early stages of an ultrarelativistic 
heavy-ion collision. 

PACS: 12.38.Aw; 11.55.-q; 11.15.Kc; 12.38.Bx 
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Typeset using REVT^ 



I. INTRODUCTION 



Recently, a lot of effort is made to study the production and equilibration of the quark- 
gluon plasma in ultrarelativistic heavy-ion collisions at RHIC and LHC 0]. Such a state 
of matter is predicted by lattice QCD calculations at high temperatures and high densities 
0. The major problems one encounters in such a study is how quarks and gluons are 
formed in these experiments and how their phase-space distribution evolves in space-time 
with collisions among the partons taken into account. At ultrarelativistic energies, the two 
nuclei are highly Lorentz contracted. When they pass through each other, a chromoelectric 
field is formed due to the exchange of soft gluons This is a natural extension of the 

color flux-tube or the string model which are widely applied to high energy pp, e'^e~ and 
pA collisions This chromoelectric field polarizes the QCD vacuum which leads to the 

production of quark/antiquark pairs and gluons via a Schwinger-like mechanism. These 
quarks and gluons collide with each other to form a thermalized quark-gluon plasma. The 
space-time evolution of the formed partons can be studied by solving the relativistic non- 
abelian transport equations for quarks and gluons As the chromofield exchanges color 

with quarks and gluons, color is a dynamic quantity. The time evolution of the classical 
color charge follows Wong's equations 



-r'u^Q'A'^ (1) 

and 

^ = Q'^F^f'^'u,. (2) 

Here, the first equation describes the precession of the color charge in the presence of 
a classical background chromofield. The second equation is the non-abelian version of the 
Lorentz-force equation. Taking Wong's equations into account, one finds the relativistic 
non-abelian transport equation PJTO[| which for quarks reads for example as: 



p,d^ + gQ'^F^ydi: + gr'V^'Ayd^Ql fix, p,Q) = C + S. (3) 

Note that there are different transport equations for quarks, antiquarks, and gluons. 
In the equations mentioned above, f{x,p,Q) is the single-particle distribution function of 
quarks, antiquarks, and gluons respectively in the 14-dimensional extended phase space of 
coordinate, momentum and color in the SU{3) gauge group. The first term on the LHS 
corresponds to common convective flow, the second term is the non-abelian generalization 
of the Lorentz force term and the third term represents the precession of the color charge 
in the presence of the background field as described by Wong's equations. On the RHS, C 
denotes the collision term which describes the collisions among the partons and 5* represents 
the source term for the production of the respective particle species. The source term S 
contains the detailed information about the creation of partons from the chromofield. It is 
defined as the probability for the production of a parton per unit time per unit volume of 
the phase space. For this reason, the development of the QGP is fundamentally governed 
by this source term S. 
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The formation and equilibration of the quark-gluon plasma within the color flux-tube 
model is studied by several authors However, in all these studies, source terms for 

parton production from a constant chromoelectric field were employed, because source terms 
for parton production from a space-time dependent chromofield were not studied in literature 
before. In fact, the chromofield acquires a strong space-time dependence as soon as it starts 
producing partons. Also the acceleration, collision and precession terms present in the 
transport equation make the field space-time dependent. This implies that the source term 
for particle production from a constant field is not applicable under these circumstances. 

We mention here that particle production from the classical field via vacuum polarization 
is studied in two different cases. In the presence of a constant uniform background field 



the particle production is computed through the Schwinger mechanism |T2[ which is an 
exact one- loop non-perturbative result. This result can also be understood as semi-classical 
tunneling across the mass gap |T^. However, for a space-time dependent field, particles can 
be produced directly by a perturbative mechanism [|T^,Q. It is observed in numerical studies 
P that the chromofield acquires a strong space-time dependence as soon as it produces 
partons. Due to these reasons, formulas based on constant fields are not applicable to such 
problems and source terms for parton production from a space-time dependent chromofield 
are necessary. In this paper, we will derive the source terms for the creation of qq and gg 
pairs from a space-time dependent chromofield in order of S**^^-*. 

The production of qq pairs from a space-time dependent non-abelian field is similar to the 
creation of e^e~ pairs from an abelian field. Apart from the color factors, the interaction of 
a quantized Dirac field with the classical potential is the same in both cases. So, the results 
obtained for the probability of creation of e^e~ pairs from a Maxwell field can be transfered 
to the creation of qq pairs from a Yang-Mills field. Contrary to that, the production of gluons 
from a classical external Yang-Mills field has no analogue in the abelian case. There exist 
direct interactions between the quantized non-abelian field and the classical non-abelian 
potential which then lead to the production of gluons from the QCD vacuum. In this paper, 
we will derive the source term for the production of gluons from a space-time dependent 
chromofield through the application of the background field method of QCD which was 
developed by DeWitt and later on extended by 't Hooft . 



As already mentioned before, in the color flux-tube model the creation of the chromofield 
is via soft gluon exchange when two nuclei cross each other 0-0]. It is the extension of the 
string model widely applied to high-energy e"'"e~-, pp- and pA- collisions P,^. However, this 
is not the only model where a chromofield is created in high energy heavy-ion collisions. 
The McLerran-Venugopalan model also predicts the existence of a chromofield at low x and 
low transverse momentum [|T^,|T^. Also recent work by Makhlin and Surdotovich \T^\ and 



many others, e.g. p0| , |21|| seems to indicate that the highly gluon dominated environment 
in heavy-ion collisions will be ruled by scales resulting from the presence of background 
fields. In any case, once there is a chromofield it will polarize the QCD vacuum and will 
produce quarks and gluons which will then form a thermalized quark-gluon plasma. In 
the McLerran-Venugopalan model the parton production is described as classical radiation 
whereas we describe parton production via vacuum polarization. It might be interesting 
to incorporate gluon creation both by classical radiation and vacuum-polarization in the 
transport equation to study the evolution of the quark-gluon plasma. 

The paper is organized as follows: In chapters II and HI we present the source terms for 
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quark and gluon production from a space-time dependent chromofield. In chapter IV, we 
analyze the production of partons from a purely time dependent background field. Finally, 
we summarize and discuss our results in chapter V and conclude in chapter VI. There, one 
can also find a treatment on the applicability of the source terms presented in this paper. 
In the end, there follows an Appendix incorporating more mathematical details than were 
suitable for the main body of the paper. 



II. SOURCE TERM FOR QUARK PRODUCTION FROM A SPACE-TIME 

DEPENDENT CHROMOFIELD 

The amplitude for the lowest order process A — > q{ki)q{k2) (see Fig. (1)) contributing 
to the production of qq pairs from a space-time dependent classical non-abelian field A""^ 
via vacuum-polarization is given by: 

M =< q{k{)q{k2)\S^^^\^ >= u\k,)iV%^A^''iK)v\k2) (4) 

with the interaction vertex: 

iV%, = igi,^. (5) 

Here ki and stand for the four momenta of the outgoing quark and antiquark respec- 
tively and A°'^{K) is the Fourier transform of the space-time dependent chromofield A°'^{x) 
with K = ki + k2: 

Al{K) = ^ / d'xA^ix) e+^^-. (6) 

and represent the Dirac spinors of the outgoing quark and the outgoing antiquark, 
respectively. Note that in the end the integrations over K have to be carried out over the 
correct kinematical region, because real quarks are to be produced. It is given by {K)"^ > 4m^ 
and > 0. This fact is not explicitly mentioned in the formulas, but has got to be kept 
in mind. 

In order to obtain the probability, the absolute square of the amplitude has got to be 
integrated over the phase space of the outgoing particles: 

W=J:[ ''\\ MM*. (7) 



spin 



Putting in the definitions for the amplitudes and performing the spin-sum, as well as 
carrying out one of the d'^k integrations yields ||14|| : 

[kK + - ^g""] (8) 

where oji 2 is defined as: 
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^1,2 = ^\ki,2? + rn^- (9) 
Using the Fourier transform Eq.(P) of the classical field in Eg.® one finds: 



From the resulting expression ([lOD, it is then possible to extract the gauge invariant 
source term for the production of qq pairs. It should be mentioned here that the probability 
of pair production Wg^^ is a real quantity because T is real (see Eq. (|^). Therefore it can be 
shown that the imaginary part of the above expression vanishes. However, for mathematical 
convenience we always keep the full expression Eq. (|10[) and take the real part at the end of 
our calculations. Again, this will not be denoted in the formulas, but must not be forgotten. 

Starting from Eq.(|iy) we find the source term for qq production: 



A^ix) J d'x2 A:{x2) J^Jd'K 5^'\K - k ~ k,) 



d^x(Pk 2(27r)6cj '7 " J 2cu2 



^iK{x-X2) ^j^^j^u ^ ^M^r. _ (11) 



The d'^K integration in Eq.(|TT]) can be carried out so that one finds: 

dWqq^ _ g"^ Aaf^\Ak-x f ji^ Aa / \ -ik-X2 f d^k2 



d^xd^k 2{2n)^uj ' J ^ ^' J 2^2 



The factors of /c2 in Eq. (|l2D can be replaced by differentiations with respect to i{xi — X2) 
which leads to: 



d^xd^k 2(27r)6^ 



^ - — ^ Al{x)e'^-'' j d'^X2 Al{x2) e-*^-^2 



[P--^— + ^ ^ F -{m' + k. ^ )gn I ^ e^'^-^---). (13) 

tO{X — X2)u 10[X — X2)iM %0[X — X2) J 2,002 

The d^k2 integration can be performed analytically which yields 

i^l = (14) 



Here Ki{z) stands for the modified Bessel function of the third kind and order one. 
Substituting Eq.(|l3D into Eq.(p!3|) and performing the differentiations Q we obtain the general 
result 



^For the differentiation of the Bessel functions see e.g. Ref. 12 
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- X2Y + {x- X2Yk'' + k-{x- X2)g'"') 
_Ko{m^-{x - X2Y)m^ -{x - X2Y + 2Ki{m^J-{x - X2)^) 



-mV- V ^ ^ ]. (15) 



-(x - X2Y 



This is the source term for the production of qq pairs from a space-time dependent 
chromofield. This source term contains all the information about the quark production from 
a space-time dependent chromofield 74"(a;) in the phase space of coordinate x and momentum 
k to the order S^^\ 



III. SOURCE TERM FOR GLUON PRODUCTION FROM A SPACE-TIME 

DEPENDENT CHROMOFIELD 

For the gluons in QCD the generating functional is given by: 

Z[J, = J [dA] [dx] [dx^] exp{zS[A, x,x^] + J ■ A + x ■ + ■ 0- (16) 

Here J, ^, and are external sources for the gauge field A and the Faddeev-Popov ghosts 
fields x^ cind x, respectively. The action S is defined as follows: 

S = J d^xC (17) 

where C consists of three terms C = Cq + Cqf + ^fp- j^g is the lagrangian density of 
the gauge field: 

£g = -\f;,F'^^'' (18) 

where is the non-abelian field tensor defined as: = 9^v4,^ — di^A^ + gf"-'"^A^^Al 
with the antisymmetric structure constant of the gauge group /'^'"^ and the coupling 
constant g. The contribution of the gauge fixing term is given by: 

^GF = -^(Gn' (19) 

where is a arbitrary form linear in the A-field that fixes the gauge. For the common 
choice = d'^A'J^ the Faddeev-Popov ghost term becomes: Cpp = ~{9fixD{.D'^x)a where 

denotes the covariant derivative = + gT°'A°'^ with the generators of the gauge 
group in adjoint representation T". The scalar products in Eq. (0) are defined as: 

J ■A = j d^xJ^A^". (20) 
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In the background field method of QCD, the gauge field is split into two parts: 

A^^A^ + Q^. (21) 

From here on, is a classical, i. e. non-quantized background field and is a quantum 
field representing the gluons. The lagrangian density of the gauge field now becomes a 
functional of A and Q: 

I^G = -\f^,AA + Q]F^'''[A + Q]. (22) 

Additionally, a special gauge is chosen. Only a gauge for Q has got to be fixed, as only 
this field is going to be quantized. The so called background field gauge is given by p6|] : 

G« = d^Q""^ + gT^^A^^'Ql = {D^[A]Q''^). (23) 

Due to this choice, the ghost term becomes: 

Cfp = -{D,[A]xi){D^[A + Q]xn- (24) 

The total lagrangian density C in the background field formalism is gauge invariant with 



respect to the type I gauge transformations WM: 



A^^A'^ = UA^U-' - '-id,U)U-\ (25) 



and 



Q,^Q'=UQ,U-\ (26) 



X^X' = UxU-\ (27) 



even after quantization. In the following, all calculations are to be performed in the 
background field Feynman gauge, i.e. with a = 1. 

For the gluon pair production from a space-time dependent chromofield in the order of 
S^^^ the vacuum-polarization amplitude is given by: 

M =< kik2\S^^^\0> . (28) 

To this order those terms of the interaction lagrangian £/ involving two Q-fields are 
given by: 

^i'^ = -lF^u[A]gr'V''Q'"' - lid.Q: - d,Ql)gr'\A'^Q'^ + Q'^A^'^) 

-\g'r''r'''{A'^Ql + Q',Al){A''^^Q^'^ + Q^'^A^''^) 

-OxQ^^gT^'A^Q'"' - ^^g^f^'f^'"' A^Q'^aIQ"'''. (29) 
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From these terms one can read off the necessary Feynman rules. For the production of 
two gluons by couphng to the field once (see Fig. (2a)), one obtains the vertex: 

{ViaTX = 9r''[-29,pK - g.pih - h), + 2g^,K;\, (30) 



where K = ki + k2. If two couplings to the field are involved (see Fig. (2b)), one receives: 

\abcd 



{V2A)tt = -w'lr'T'ia.xgup - g,pgux + g,ugxp) 



+r'^T^\gpugxp - gpxgup - g^^pgux) 

+r''r'\gpugxp-g,pgux)]. (31) 



The above rules coincide with those derived in |^ . The total amplitude for the produc- 
tion of gluon pairs is now equal to the sum of the single contributions: 

M = MiA + M2A- (32) 

If the correct weight factors are included, so that the lagrangian density can again be 
retrieved, one finds for the first term: 

M,A=^^J d'K6^'\K -h- k2)A^^{K)e'^{h)e'p{k2){V,A)t.% (33) 

where finally the allowed kinematical region of integration is given by {K)"^ > and 
> 0. The second term is given by: 

M2A = \j d%d%6^^\k, + k2-ks~ k,)A''^ik,)A'^\k,)e"'iki)e''P{k2)iV2ArX^ (34) 

where the limitations {k^ + k^)^ > and k^ + k^ > have to be obeyed in the end. The 
probability can again be obtained by making use of Eq.(|^. In order to obtain the physical 
polarization of the gluons, we use the following spin-sums: 

E e^'{h)e*^''ik^) = e^(A:2)e*^'(A:2) = -g"''' , (35) 

spin spin 

but later on deduct the corresponding probabilities for the processes involving ghost 
instead of gluons as shown in Fig. (3). In our case, the total probability for the production 
of gluons not yet corrected for ghosts consists of four terms which add up: 

= WiA,lA + WiA,2A + W2A,IA + W2A,2A- (36) 

By definition, we have for the first term in Eq. 



1 /• d^ki d^k2 



For the integration over the phase-space of the gluons there is an additional factor of 
1/2 because the particles in the final state are identical. The introduction of the Fourier 
transform Eq.(^) and the application of 
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j^abd j^a'bd g^aa' (38) 

yields: 

K27r)6 J 1^14 

[2gf,f,'{h + k2)^ - 2iki + k2)^{ki + k^)^' + {k^ - A;2)m(^i - ^2)/.']- (39) 

Again, it has to be noted that the real part of this quantity has to be taken. By the 
same means as in the quark case, we obtain the first contribution to the gluon source term: 



d^xcPk 8(27r)6 7 k^k^ 
[2(7^^' (A; + k2f - 2{k + k2)^{k + k2)^' + {k- k2)^{k - ^2)^']. (40) 

Making use of the following relation 

<^ ik2-{x-x') ^ 8^ (4^1) 

Ihl (x-x')2' ^ ^ 

and explicitly calculating the derivatives afterwards results in the final form (see Ap- 
pendix) : 

^ {x — x'y ^ (x — x')"^ {x — x'y (x — x')^ {x — x'y 

Now we perform the same procedure with the other contributions. First for the two 
mixed contributions we get: 

1 r d^k d^k 
= W..,,, = - / ^^M..Mr,. (43) 

After expressing the probability in terms of the Fourier transforms of the A-fields: 

WiA,2A = W2A,1A = 

-ig^ f d^ki d^k2 



I ^^d%d^k4d^Kd^X3d^x^d''x6'^^\K - k^ - k4)6^^\K - h - /ts) 

J rh-\ /Co 



64(27r)io J kl kl 



we can extract the source term: 
W = 8(Sf / '^/kAd^Kd^x,d^x,8^HK - k, - k,)5i^\K - - fc^) 

At this point, we again need Eq. (^l]) to find (see Appendix): 



dW2AM _ 3i^' f d':,'e^H^'-^)A^\x'){A'^{x') ■ A^' {x))r'^%—^ + 2^^^^^). (46) 



And finally for the last term: From the definition 

«W.4i;/p||^p|i^M..M.-. (47) 

we obtain 

W,A,2A = ^^1^ / '^'^d%d'kd%d%d'x,d'xd%d'x',e^^'--^^^^^^^ 

6^^\ki + k2-ks- ki)5^^\k^ + ki-k'^- k'^)A''^'{x^)A^^{xi)A'''^'\x'^)A'''^\x'^) 

[O^ „ ( fabx fxcd I radx rxcb\ / ra'bx' fx'c' d i ra'dx' rx'c'b\ i O/l ^ ^ racx ra'c'xl / ac 
i^fl'AiAVA'U / + J J >[J J +J J ) + ^^9fiii'gx\'J J J- 

And from there the source term: 

32(27r)i4i PA;^ 3434434 
5^^\k + k2-k3- ki)5^^\k^ + ki-k'^- k'^)A''^'{x)A''^{xi)A'''^''{x'^)A''^\x'^) 

^ ( fabx fxcd I fadxfxcb\ffa'bx'fx'c'd i ra'dx' fx'c'b\ i 10^ ^ racx ra'c'xl t Ar\\ 

[gtj.xgfi'x'[j J + J J )[j J +j J > + i'^gfj.^i'gxx'j j \- i4yj 

Here we apply Eq. (^Tj) to receive (see Appendix) 

r„ „ I eabx fxcd , radx fxcb\ / ca'bx' fx' c'd , ra'dx' fx'c'b\ , io„ „ racx ra'c'x} 

[fl'AiAVA'u J + J J >[j J +j J ) + i-^gfiti'gxx'j J \- [ou) 

As already remarked before, we also have to calculate the corresponding terms for the 
ghosts in order to be able to remove the contributions by unphysical polarizations of the 
gluons. To the ghost matrix element, there contribute the two Feynman diagrams shown in 
Fig. (3) with the interaction vertices 



and 



ivf/rj" = +gr'\k, - k^),, (51) 



(v/aV = -^g'g,xir'^r' + r'T"). (52) 



These vertices are obtained from the ghost lagrangian density: 

= -id,x''^)r'"'Kx'' - {r'"'Ai[A]x'^)d,x^ - {r'''4x'^)r''^A'^x' (53) 

Including the correct weight factors we obtain the amplitude: 

^M^pyd ^ (Mf/)^'^ + (Mf/)'*^ (54) 

with 



10 



and 

{M^lf' = \j d%d%6^^\h + k2-h- h)A'''{h)A'\h){V,Y)X- (56) 

The probability is defined as for tlie gluons but without the factor 1/2 for the phase 
space: 

^ = y (2^0^ (2^^^^ ) ) . (57) 

This time, there are only two contributions, because the cross-terms vanish: 

W = WZa + WZa- (58) 
Starting from the definition of the first term 

tFP f d^ki d^k2 /■ j,jFP\bd/- i\jFP\*bd 



\2 



we carry out the same steps as for the derivation of Eq. (^) from Eq. (p7D through 
Eq.®. From 



Q„2 r d^k d^k 



we extract the source term: 

dW,Z^ 3g' 



d^xd^k 4(27r)5A;0 



Analogously, we obtain the second contribution from its definition 

d^ki d?k2 
(2^0W(2^032A| 



^2%A = I Jlto/.to iMiinMi/r^ (62) 



and finally get: 

^2%A = / ^^d%d%d%d%5^'\k, + k2-ks- k,)5^^\h + k2-k',- k',) 

A^^'{k^)A^^{ki)A*'''^'\k'^)A''''^'{k'^)g^Xg^,'X'{^^T'^ + f^dx jxcb^^ja'bx' jx'c'd ^ ^a'dx' ^x'c'by ^g3^ 

Comparison to Eq. (pDD yields the source term: 



^ ^ ( fabx fxcd I radx rxbc\ / ra bx rx c d 1 ea dx fx be \ f(iA\ 

gtj.xgtM'x'[j J + J J )[j J +j J >■ it>4j 
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In the end, we obtain the required source term for the production of gluon pairs from a 
space-time dependent chromofield via vacuum-polarization by combining all the contribu- 
tions in the following manner: 



Wgg = W,AM - W^lXlA + 2W^2A,1A + W2A,2A ' W,%a- (65) 

So, we receive: 



^^^^ - ^ /rf^xV'=-(^--')^-^{^/A«'^(x)A«'^'(xO[3fc,V-8^,,,rz^'^ ""'^ 



'\2 



d^xd^k (27r)5A;0y {x-x'^'A" v ^ mm ^mm (^^ _ 

(x-x')^' , , .{x-x')f, g^f,' (x - x')f,{x - x')i,' 



+5( Vt + ji) + 6H^^ - 12 



X — x'Y (x — x'Y — x'Y ~ ^'Y 

-?,ig^A''^{x')A'\x')A'^'^^'{x)r'-'^K^g,,, 
-^c/^A"'^(x)A"^(x)A"'^'(x')A"'^'(x') 

This is the source term for the production of gluon pairs for any arbitrary space-time 
dependent chromofield in the order 5**^^^ 

It can be checked that this source term is gauge invariant with respect to type-(I)-gauge 
transformations in the following manner (see also |2^). One can observe that the gauge 
invariant part of the lagrangian which in each term contains two Q fields or two ghost fields 
is given by: 

C = -\gr'' F^^MWQ'" 

-^[DMW? - {DM\xi)DM\x''. (67) 

This lagrangian density is gauge invariant with respect to the gauge transformations in 
Eqs.([25|), (^), and (^). It can be checked that the free terms, i.e. terms not involving 
A, of the gauge invariant lagrangian density (Eq. (|67D), when evaluated between an initial 
vacuum state |0 > and a final physical two gluon state < ki^k2\, do not contribute to the 
probability. Hence, the expression we have used for C in Eq. ( [29| ) together with Eq. (^) for 
the ghosts gives the correct result for the gluon pair production probability, invariant with 
respect to the above gauge transformations. 



IV. SOURCE TERM FOR PARTON PRODUCTION FROM A PURELY TIME 

DEPENDENT CHROMOFIELD 

Our main purpose is to use the obtained source terms of quarks and gluons in the 
relativistic non-abelian transport equations to study the production and equilibration of the 
quark-gluon plasma expected to be formed in ultra relativistic heavy- ion collisions (URHIC). 
Solving the transport equations with all the effects, such as collisions among the partons. 
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acceleration of the partons by the background field, precession of the color charge in the 
group space and production of the partons from the background field involves much more 
numerical work as was done in . In the following we are going to analyze an example: We 
will consider a time dependent chromofield and examine the behavior of the source term in 
the phase-space. Solutions of the transport equations including these source terms to study 
the production and equilibration of quark-gluon plasma will be presented elsewhere. 



A. Source Terms for any Arbitrary Time Dependent Chromofield 

For a purely time-dependent field Eq. (pJ]) reduces to: 
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J dt, Al{t^) J^Jd'K 5^'\K -k-k,) 6^'\K) 



d^xd^k 2(27r)3cu 'J ' "'^'J 2uj2 

After integrating over d^K and d^k2 and taking care of the Dirac 6 factors we find 
dW - r p2iaj{t-t2) 



dt2^ [2Lu'{A''{t) ■ A''{t2)) - 2{A''{t) ■ A:)(A»(t2) ■ k)]. (69) 



d^xd^k 2{2'KfujJ 2uj 

Here, we have presented Eq. (|69D in the form of scalar products of three-vectors. The 
last remaining time integration then constitutes a Fourier transform from t2 to 2uji. The 
final result for the source term in the case of a time-dependent chromofield is 

dWqq _ g^^/2^ ^2iuU,2(7a . . 



d^xd^k 2(27r)3u;2 



uj\A\t) ■ A*\2uj)) - {A\t) ■ k){A*''{2uo) ■ k)]. (70) 



Analogously, the Eqs. (0,|5, and ^) for the gluons and the Eqs. (^ and for the 
ghosts simplify to (see Appendix): 



dWiAM _ 3gW2n 



^Itk^t [(^a(^) . - )(^*a(2fcO) . - ) _ 2[A\t) ■ A*\2k^))l (71) 



d^xd^k 2(27r)3 ^ ' k"^'^ ' ' 



dW2A,iA -3ig 



3 



d^xd^k 4(27r)3A;0 
dW2A,2A g"^ 



I dhe^'^'^''-'^A-^{t,)A'\h)A'^'^\t)r'^'^g^^,, (72) 



d^xd^k 32(27r)3(A;0)2 

„ ( fabx fxcd I fadxfxcb\/fa'bx'rx'c'd i ra'dx' rx'c'b\ i 10^ ^ facx ra'c'xl /''7Q^ 

[g^J.xg^l'X'[J J + J J ){J J +J J ) + ^2g^^igx\iJ J J, {(A) 
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and 



^ „ ( fabx rxcd < radx rxcb\ ( ea'bx' ex' c' d < ra'dx' fx'c'b\ ('7r.\ 

gi,\9t,'\'[j J + J J ){J J +J J )■ (75J 
Adding the gluon and deducting the ghost contributions, we obtain: 



dW,, _ ^2..o,^ 3l ^^^^^^^ _ A)(1-(2A;°) . A) _ 4(l'^(t) ■ i-(2A:°))] 



+ * A^°)*(2A;°)A'^'^'(t)r'-^?^^, 

^ ( fabx cxcd < radx fxcb\ / ra'bx' rx'c' d i ra'dx' fx' c'b\ i o/l ^ ^ facx ra'c'xl^ (^Ci\ 

[gfi\gfi'x'[j J + J J )[j J +j J > + ^^gfifi'gxx'j j Jj- Uoj 

This is the general source term for the production of gluon pairs in the presence of a 
purely time dependent chromofield. 



B. A Special Case 

We now choose a special form of the field in order to get some insight into the behavior 
of the obtained source terms for quarks and gluons. For simplicity we choose the field to be 

A'^^it) = A,„e-l*l/*», to > 0, a = 1, 8, (77) 

and all other components are equal to zero. Many other forms could have been taken. 
We have chosen this option just to get a feeling for how the source term in the phase-space 
behaves. This choice is also inspired from the numerical study which shows that the 



decay of the field is close to this behavior. In any case the actual form of the decay of 
the classical field can only be determined from a self consistent solution of the relativistic 
non-abelian transport equations, as the transport equations have to be combined with the 
energy momentum conservation equations for particles and fields. That is why the exact 
form for the decay of the fields (due to the production, acceleration of partons, precession 
of color, collision among partons, and expansion of the system) can be determined from 
self consistent transport studies (see [2^] for details). We will use the general source terms 



(see section III and IV) in the relativistic transport equations to describe the production 
and equilibration of the quark-gluon plasma in ultra relativistic heavy-ion collisions in the 
future. 

The Fourier transform of Eq. ([77| ) is given by 



A"3(-2|A;|) = A,„^ ^-^ . (78) 

Putting Eqs. (^) and (|78D into Eq. (^) and afterwards summing over color space yields: 
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d^xd^k (27r)2 ' ™' 1 + 4cj2t2 ^2 ' 

with m|i = + and where fcr is the transverse momentum. Using, in the same way, 
Eqs. and (|7|) in Eq.@ results in: 



^^gg ^ 24as n2 2tf't -|t|/to / o _ - 



+ 27r ^^"^ l + (A;0)2t2(fcO)2- 

The contribution of the mixed term (HBI) vanishes for all fields of the form A"-^{x) = 
^infi^)- For our choice the contribution of the last term in Eq. (|50D is zero. 

In the next chapter, we will present plots for the source terms derived above (Eqs. (|79D 
and (|80|)) for our choice of input parameters. 



V. RESULTS AND DISCUSSION 

The general results for the source terms for the production of quark-antiquark pairs 
and gluon pairs by vacuum polarization in the presence of a general classical space-time 
dependent chromofield to the order S^^^ are given by Eq . ([15|) and Eq . (|66|) , respectively. 
Continuing from this point, i. e. using these source terms in relativistic non-abelian transport 
equations to study the production and equilibration of the quark-gluon plasma, involves 



extensive numerical computations |2J] which have to be represented later on. In order to 
get an idea about the behavior of the source terms, we derived the general expressions for 
them in the presence of a purely time-dependent field. Subsequently, we assumed a certain, 
physically motivated form for the background field, given in Eq. ([77| ) to obtain Eqs. ( |79[ ) 
and (|80|) . We now present the results from the above equations. 



We relate the energy of the produced particles to the rapidity y by: 

k° = krcoshiy), (81) 

where kr is the transverse momentum. If not stated otherwise, we choose the following 
parameters: as = 0.15, Aj„ = 1.5GeV, kx = l.5GeV, y = 0, and to = 0.5/m and the quarks 
are considered to be massless. These values might not correspond to an exact combination 
of values in an URHIC. We have just chosen a set of values as an example in order to 
demonstrate the properties of our source term. 

In Fig. (4), we plot the source term 

9- (82) 
^ - Wfc- ^^^^ 

for quarks (Eq. (^)) and for gluons (Eq. (pO[)) versus time for the above choice of 
parameters. All quark graphs are to be multiplied by a factor of two. The exponential 
decay originates mostly from the decay of the model-field. Of course, every decaying field 
will transfer this quality to the production rate, but this is not necessarily an exponential 
decrease. The oscillatory behavior is due to the exponential factor with imaginary exponent. 
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This factor is already present in the general formula. Physically, it seems to indicate that 
there exist periods of particle creation and particle annihilation which follow each other 
periodically. Additionally, the frequency of this oscillation is varying with the energy of the 
produced particles. This oscillatory behavior of the source term will play a crucial role once 
it is included in a self consistent transport calculation. It can also be seen in the Figure that 
there are considerably more gluons produced than quarks. 

In order to get a better overall understanding, we present in the following two Figs. ((5) 
and (6)) the source term integrated over positive times, i.e: 

- (83) 



2d?xd?k 



The time-integrated source terms can be regarded as a measure for the net-production 
of particles in an infinitesimal volume around any given point in the phase-space. It does 
not show the oscillatory behavior which gives a totally different picture for different times. 
For our choice of the field (Eq. (|77D ) T takes the form: 

-^^i^M.n)\ ^ , !: ..j ^^Vt (84) 



and 

dWgg _o «5 ^0 



2(Pxd?k (27r)4' ™' ' ^l + A{kTYcosh^{y)tl' ' cosh^^/) 

^^^^^'^^^ {kT? cosh^(y) ^ 1 + {krY lo^h\y)e ^^^^ 



respectively. 

The decay behavior with the transverse momentum observed in Fig. (5) is mostly due 
to the choice of the field. Only in the second contribution to the gluon source term there is 
already a factor l/{k^)'^ present in the general formula. Choosing e.g. a Gaussian form for 
the field, would also have resulted in a Gaussian decrease of the production rate with the 
transverse momentum of the particles. As the momentum structure of the general equations 
is mostly based on the fc°-component, the origin for the typical rapidity behavior is mainly 
the same as for the behavior for changing transverse momentum, see Fig. (6). 

In the last plot Fig. (7), we plot the ratio of the source terms for quark-antiquark pairs 
and gluon pairs: 

^ dW,q dWgg 

2(Px(Pk' (Px(Pk' ^ ' 

It can now be seen directly from Fig. (7) that there are less quarks produced than gluons. 
We observe for this model field that a stronger coupling, a stronger chromofield and/or a 
slowlier varying field emphasize the gluon-pair production even more in comparison to the 
production of qq pairs. 
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VI. CONCLUSION 



We have derived the source term for quark- ant iquark and gluon pair production via 
vacuum polarization in the presence of a space-time dependent chromofield in the order 
S^^\ We have used the background field method of QCD for this purpose. The production 
of quark-antiquark pairs from a Yang-Mills field is very similar to the production of electron- 
positron pairs from a Maxwell field in QED. To the order S^^^ the source term for gluon pairs 
in the presence of a general space-time dependent chromofield consists of three contributions, 
each of a different order in the coupling constant g. 

To obtain an insight into the source terms, we have derived them for a purely time 
dependent chromofield and have presented the results for an exponentially decaying field. 
In comparison to the arbitrary space-time dependent field, the equations obtained in this 
case are considerably simpler but still too complicated as to interpret them at first sight. 
For this reason, not for giving a realistic physical prediction, we evaluated the source term 
in the purely time-dependent case for a given form of the chromofield. It was chosen, so 
that analytical solutions were possible, but still it is motivated by physical arguments. The 
results thus obtained have been plotted and discussed. 

We observe that the source term not only creates partons from the field but also destroys 
them to produce field. One can recognize periodic fluctuations of the production rate with 
time which become more rapid for higher particle energies. This behavior is inherent to the 
general time-dependent formulas and is not connected to a special choice of the chromofield. 
This feature may play an important role in the evolution of the quark-gluon plasma. The 
mere existence of the field might complicate the transport studies. For example, one usually 
solves the hydrodynamic evolution equation |26] when the plasma is in local thermal equilib- 
rium. The presence of the chromofield results in the necessity to solve chromohydrodynamic 
evolution equations, thus one needs to know how fast the field decays. The description of 
all processes will become substantially easier if the field strength is already negligible before 
equilibration sets in. In all plots it can be seen that there are considerably more gluon pairs 
than quark-antiquark pairs produced by this mechanism. 

Finally, we discuss the range of validity of the source terms derived in this paper. Being 
perturbative, our result is only applicable in situations, where certain conditions are satisfied. 
Firstly, in the framework of our approach we can merely properly describe the production 
of partons with momentum p larger than the field strength A multiplied by the coupling 
constant g. This means, if the field is very strong and gA is truely large, our calculation is not 
valid, even if the coupling g is weak. This is similarly true for the Schwinger mechanism for 
pair production in QED in the presence of a weak field (see Eq. 6.33 in [|l^]). Although QED 
is a weak coupling theory the perturbative calculation of Schwinger is applicable only in the 
limit of small eA. However, for a very small coupling, a truely large background field has 
got to be present before this condition is violated for reasonable momenta p. Secondly, for 
the perturbative method to be applicable gA has got to be larger than the non-perturbative 
scale Aqcd- For the condition gA < Aqcd the perturbative calculation is not valid and 
parton production has to be computed via non-perturbative methods, which is unfortunately 
non-trivial. So, the region in which our perturbative approach is applicable is given by : 
p> gA> Aqcd- 

Let us consider a realistic situation in a heavy-ion collision at LHC (Pb-Pb at ^/s= 5.5 
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TeV). For the energy density formed in this experiment one expects e ~ 1000 GeV/fm' 



3 



|, |28|jl8| and for the strong couphng constant as ~ 0.15 |23]. Assuming that all the energy 
density e is deposited in the field sector we can get a rough estimate for gA {E ~ ^/2e and 
gA^ ~ E) to be: gA ~ 2 GeV. Hence, initially, a perturbative treatment might be applicable 
for the production of partons whose momentum is greater than 2 GeV or so at LHC. At 
RHIC (Au-Au collisions at ^/s = 200 GeV) the typical energy density and coupling contstant 
would be e ~ 50 GeV/fm^ and as ~ 0.33 respectively. For this situation one would expect 
roughly gA ~ 1 GeV. So, for these input parameters at RHIC and LHC, the minimum 
momentum above which one would apply perturbative calculations is about 1 and 2 GeV 
respectively. It is interesting to note that these values are more or less equal to the minium 



momentum cut-offs used to compute minijets by using pQCD at RHIC and LHC |l29| , |28 



We have to mention that even if gA is small the contributions of both diagrams have to be 
taken in order to maintain exact gauge invariance. 

It is clear from the above arguments that during the very early stage of heavy-ion col- 
lisions our perturbative calculation is applicable for the production of partons having a 
momentum greater than gA. Begining from such values of the field strength, it is now nec- 
essary to investigate how the field decays owing to all effects such as particle production 
from the field, acceleration of the partons by the field, precession of the color charge in the 
presence of the field, collision among the partons, and expansion of the system. As time 
progresses the field strength decreases due to the effects mentioned above. However, as long 
as gA is greater than Aqcd, our perturbative method is applicable. The decrease of the 
product gA with time allows for the description of partons with less minimum momentum 
than right at the start of the URHIC. The picture breaks down however, when gA reaches 
Aqcd- In that case non-perturbative calculations for parton production are needed, which 
is beyond the scope of this paper. Unfortunately, when trying to perform non-perturbative 
calculations one always encounters many more obstacles. 

How rapidly the boundary of applicability of our picture is reached and how the quark- 
gluon plasma evolves (up to a time where the condition gA > Aqcd is violated) with 
our treatment of parton production from the space-time dependent chromofield, must - 
preferrably - be determined from a self-consisitent transport calculation with all the above 
effects taken into account. This project will involve a considerable amount of numerical 
work and will be undertaken in future. 
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APPENDIX 



Gluon Source Term 



We start from definition (^) and introduce the expression (|33|) and subsequently Eq.(pOD 
to obtain: 

f-bdfa'M^_2g^^K, + Qupih - k^), + 2g^,Kp][-2g^,pKl + g,p{k[ - k'^)^, + 2^.^;] (87) 

Now we make use of tlie Fourier transform of Eq.(|^). In the next steps, we first carry 
out the integration over d'^K' and then over d'^K: 

[2^MM'(^i + ^2)^ - 2{ki + k2)f,{ki + /C2)^' + {ki - k2)^{ki - k2)^'] (88) 



In the previous step, use has been made of Eq.(p8D. Now, we extract the expression (^0]) 
for the source term. From there, we proceed by replacing the occurrences of /c2 by derivatives 
with respect to i{x — x'): 

dW^AM fd^:,'^e^H^-^U'^^{x)A'^^'{x')[2g,,ik+ ^ 



d^xd^k 8{27i)^k^ J kl V n^MMv ' id(x - x')' 

-2{k + ^ ),ik + ^ + {k- ^ ),{k - ^ ,J ,.]e-^^-(--') (89) 

io[x — x') io[x — x') io[x — x') io[x — x') 

Now the integration over d^k2 can be carried out with the help of Eq. (^T]) In the end, 
we explicitly calculate the derivatives and obtain Eq.(|i2|). 

Now we carry out the same procedure for the other contributions. First, we calculate the 
probability of the mixed term. Starting from the definition Eq.(P3|) we use the expressions 
(0) and (|34D and get by subsequent introduction of Eqs. ( pOD and (|31|): 



W,A,2A = W2AM = / ^^d%d%d''K5^'\K -k,- k,)5^'\K -k,- k2) 

A''^'{k:,)A'\k^)A*'''^'\K)[-2g^,pK, + g,p{kr - ^2)^' + 2g,,,K^] 

+r'''r'^r'%9,u9x, - 9,x9.p - 9,p9ux) 

+r'''r''r'\9,.9xp - 9,p9ux)]- m 

Now, we contract all possible color and Lorentz indices and make use of the Fourier 
transform (^. From there the source term (|45|) can be extracted. In the next step, the d^K 
and the d'^k^ integration are carried out, followed by the d^fcs integration: 
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d^xd^k 8(27r)io7 k^k^ 

A''^{x3)A'\x^)A-'^^'{x)r''''{k + k2)xg,,> 



-Zig^ f d^k2 4 4 



8(27r)ioy k^k^ 

A"^(x3)A^'(x4)A"'^'(x)r'"^(A; + A;2)a^?;.;.'. (91) 

The new Dirac 5 distribution can be used to perform one of the integrations over d^x: 



dHd^k 8(27r)6 7 k^k^ 

Now, k2 has got to be replaced by a derivative with respect to 2(^3 — x) in order to enable 
the use of Eq. (^Tj) as a menas to obtain Eq.(^). 

The source term of the last contribution remains to be calculated. We receive by intro- 
ducing Eq. (|1]) into Eq.(||) and all that into Eq. (p|): 

^^^'^^ = 12^ / '^-^'^-^d%d'k,d%d% 
d^'^^k, + k2-ks- k^)6^'^\k3 + k^-k'^- A;;)A'^'^(fc3)A'=^(A;4)A*'^'^'(A;^)A*^'^'(A;;) 

+r''T^\g^.ugxp - m^a)] X [r'^'''r'^''^{g^,'ygup - g^'pg^y + v^^av) 
+/"''"7"''^'(v^^?A'p - gp'ygup - gp'pgux') + r''""' r''\gp'ugx'p - gp'pg.y)] (93) 

Now, all possible indices have to be contracted: 

W2A,2A = ^^1^ / '^'^,^'^'^^ d'k,d'k,d%d%d'x,d'x,d%d'x^^^^^ 

5^''\ki + k2-k3- k^)5^^\k3 + ki-k'^- k'^)A'"'{x3)A'\xi)A'''^''{x':i)A''^'{x'^) 
[^gpxgp'x'ir'^r' + + fa'dx' jx'c'b^ ^ 24g^,,gxx'r''r'''n (94) 

From this expression we get the source term (^). In the following steps the integrations 
over all d'^ks and afterwards over d'^x^ and d'^ performed: 

dW2A,2A _ g^ f d^k2 -J4, .4, / .4 .4 / ,4 > i\k:yx+(k+k2~h,)-Xi-K-x'.-(k+k2-K)-x\] 

d^xd^k ~ 32(27r)i4 J k^kf ^'"^ ^'"^ ""'"^ ""'"^ "^^^ 

(x) A'^ (X4) A'''^'' {x'^)A''^' (xl) 

r„ „ / rabx rxcd , cadx fxcb\ / ra'bx' rx' c' d , ra'dx' rx'c'b\ , -|o„ „ racx ra 

[gpxgfi'x'U J + J J )[j J +j J ) + ^2gp.p'g\\'j j 

f d^k2 



32(27r)i4 J p|^'^4rfMrfMe^^'^+'^^-^^^-^^n2vr)^5(^)(a: - x,){27rr6^^\x', - x',) 

A"'^(X)A"^(X4)A"'^'(X[5)A^'^'(X;) 
r „ / fabx fxcd , radx fxcb\ / ra'bx' rx' c'd , ra'dx' rx'c'b\ , io„ „ racxra'c'x} 

[gpxgfi'y[I J + J J )[J J +J J ) + i-2gp^igxyj J \ 

§dHe'^''+'''^'^"""='^^"''(x)A"^(x)A"'^'(x;5)A'='^'(x3) 



g 



4 <- d^k 



32(27r)6 7 fcofco 

[gpxgp'yir'T'' + f-'^-r^){f-'^-'f-'-''^ + f'i-'fc'b^ + I2g,p>gxyr''r'''''] (95) 
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For the integration over d^k2 we again make use of Eq. (^T]) and obtain the final expres- 
sion Eq.®. 



Gluon source term for the time-dependent case 

In the case were the A fields are only time dependent, we receive from Eq. (|40|) by 
integrating over the spatial coordinates: 



dW,AAA 3(?2 r d'k 



d^xd^k 8(27r)3i k^k^ 

[2g^^' {k + k2f - 2{k + k2)^{k + k^)^' + {k- k2)^{k - k2)^'] (96) 

By eliminating all Dirac 6 distributions one obtains: 



d^xd^k 8(27r)34(A;0)2 
4[-2(l'^(t) ■ l"(t'))(/t°)' + (A'^it) ■ k){A''{t') ■ fc)], (97) 

where the result (^) can be obtained incorporating Fourier transforms. 
As before, we integrate in Eq. ( ^5D over all spatial coordinates and subsequently integrate 
over all d'^ks followed by dt^: 



dW2A,iA _ -3ig'^ f d^k2 
~¥xdFk ~ 8(27r) 



-3i^3 [ d^k2 



I T:^d%d^kd^Kd%d%5^''\K -ks- ki)5^^\K -k-k2) 

J rCt rh<2 

6^^\k3)6^^\h)e'^''-^'''+^^^^^ 



8(27r)4y k^k^ 

A'^^(t3)A^"(t4)A'^''^'(t)r'"^(A: + fc2)A^7MM' 

= ^ / ^^d%d%5('\Uk2){2mt,-t,)e^^('^^^^^^^^ 

A'''^{ts)A'\U)A''''^'{t)r''''^ik + k2)xg,,, 



]{2n)%koy 
-3zg^ 



dhe^'^" -*) (ts) {t3)A^'^'' (t ) r'^^g^^' . (98) 



4(27r)3fc0 

Starting from Eq. (^9]), we carry out similar steps for the last gluon term: 

d'^xd^k ~ 32(27r)5 J k^k^ ^ K^aua^at^e 

5^'\h)S^'\k',)S^'\k',) 

6^^\k + k2-k3- k4)S^^\k3 + k^-k'^- k'^)A''^'{t)A^^(ti)A'''^'\QA'''^\Q 

r„ „ / rabx rxcd , cadx fxcb\ / ra'bx' rx' c' d , ca'dx' fx'c'b\ , io„ „ cacx ra'c'x} 

[gfj.\gfi'\'[j J + J J >[J J +J J ) + i^2gf,f,>gxx'J J \ 
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32(27r)5y k^k^ 3434 

S^'\k + k2){27r)Sit - t4)(27r)<5(4 - t'^)A''''it)A'\u)A'''^''{t',)A''^'{Q 
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Fig. 1 Feynman diagram for tlie production of an quark-antiquark pair 
by coupling to the field A once. 
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Fig. 2 Feynman diagrams for the production of two gluons by coupling 
to the field A once (a) or twice (b). 




Fig. 3 Feynman diagrams for the ghosts, corresponding to the gluon 
vertices in Fig. (2). 
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Fig. 5 The dimensionless time-integrated source terms T for the pro- 
duction of gluon pairs (sohd line) and quark- antiquark pairs (dashed 
line) versus transverse momentum kx in MeV . 
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Fig. 6 The dimensionless time-integrated source terms T for the pro- 
duction of gluon pairs (sohd line) and quark-antiquark pairs (dashed 
line) versus rapidity y. 
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Fig. 7 Ratio R of the time-integrated source terms for the production of 
gluon pairs and quark-antiquark pairs versus transverse momentum 
in MeV. From top to bottom we have first the graph for the standard 
parameters: coupling constant as = 0.15, initial field strength Ain = 
1.5GeV, and decay time to = 0.5/m. In the subsequent graphs we 
increased the coupling constant to as = 0.3, the initial field strength to 
Ain = 3GeV, and the decay time to to = l/"^> respectively. 
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